defaultEncoding <- "UTF8"
library(dplyr)
library(sf)
library(terra)
library(raster)
library(DT)
library(ggplot2)
library(plotly)
library(leaflet)
casos <-
st_read(
"/vsicurl/https://marcelocubero.github.io/capas_proyecto/casos.geojson",
quiet = TRUE
)
casos <-
casos %>%
st_transform(4326)
casos_ <-
read.csv( file= "https://marcelocubero.github.io/capas_proyecto/casos_.csv",
sep = ",", encoding = "UTF-8")
fallecidos <- rast("/vsicurl/https://marcelocubero.github.io/capas_proyecto/fallecidos.tif")
fallecidos_rl <- raster::raster(fallecidos)
leaflet() %>%
addTiles(group="OSM") %>%
addProviderTiles(providers$Esri.NatGeoWorldMap , group="NatGeo") %>%
addProviderTiles(providers$CartoDB.DarkMatter, group="CartoDB-Black") %>%
addRasterImage(
fallecidos_rl,
opacity = 0.7,
group= "Fallecidos",
) %>%
addPolygons(
data = casos,
color = "black",
fillColor = "transparent",
stroke = TRUE,
weight = 1.0,
group= "Casos"
) %>%
addLayersControl(
baseGroups = c("OSM", "NatGeo", "CartoDB-Black"),
overlayGroups = c("Fallecidos", "Casos"),
options = layersControlOptions(collapsed = FALSE)
)